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ABSTRACT 



We present new Near-Infrared {J ^K) magnitudes for 114 RR Lyrae stars in 
the globular cluster uj Gen (NGG 5139) which we combine with data from the 
literature to construct a sample of 180 RR Lyrae stars with J and K mean 
magnitudes on a common photometric system. This is presently the largest such 
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sample in any stellar system. We also present updated predictions for J,K- 
band Period- Luminosity relations for both fundamental and first-overtone RR 
Lyrae stars, based on synthetic horizontal branch models with metal abundance 
ranging from Z=0.0001 to Z=0.004. By adopting for the uj Cen variables with 
measured metal abundances an a-element enhancement of a factor of 3 (~ 0.5 
dex) with respect to iron we find a true distance modulus /xq = 13.70±0.06±0.06 
(random and systematic errors respectively), corresponding to a distance d — 
5.5 ± 0.03 ± 0.03 Kpc. Our estimate is in excellent agreement with the distance 
inferred for the echpsing binary OGLEGC-17, but differ significantly from the 
recent distance estimates based on cluster dynamics and on high amplitude 5 
Scuti stars. 

Subject headings: globular clusters: uo Centauri - stars: evolution - stars: hori- 
zontal branch - stars: oscillations - stars: variables: other 



1. Introduction 

The Period-Luminosity {PL) relation for Classical Ccpheids and, in the near-infrared, 
for RR Lyrae stars is one of the most important methods to estimate distances and to 
calibrate secondary distance indicators (Bono ct al. 2001; Freedman et al. 2001; Saha ct al. 
2001; Cassisi et al. 2004, hereinafter C04; Catelan, Pritzl, & Smith 2004, hereinafter CPS04; 
Gieren et al. 2005). Current trigonometric parallaxes from Hipparcos (Feast & Catchpole 
1997), HST (Benedict et al. 2002), and from ground based telescopes (Lane et al. 2000; 
Kervella et al. 2004) have provided a unique possibility to cahbrate the zero-point of these 
relations. Unfortunately, the number of Galactic calibrators with accurate trigonometric 
parallaxes is very limited, and for RR Lyrae stars it is restricted to a single object, namely 
RR Lyr itself. Therefore, we still lack a sound knowledge of the slope of the different PL 
relations based on purely geometrical methods. To overcome this problem, the astronomical 
community devoted a paramount observational effort to measure the pulsation properties of 
these robust distance indicators. The ideal laboratories for these studies are the Globular 
Clusters (GCs) and the dwarf spheroidals (dSphs) in the Local Group as the RR Lyrae stars 
in a given system are located at largely the same distance. 

It has been demonstrated by Longmore et al. (1990) that cluster RR Lyrae follow a 
well-defined K-band PL relation but we still lack a comprehensive empirical investigation of 
its properties. For this reason, during the last few years we undertook a detailed pulsation 
and evolutionary investigation of RR Lyrae stars (Bono et al. 2003 and references therein) 
and, at the same time, we also started a large observational project aimed at collecting new 
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Near-Infrared (NIR) data for RR Lyrae stars in Galactic (M92: Del Principe et al. 2005) 
and Large Magellanic Cloud clusters (Reticulum: Dall'Ora et al. 2004). 

Although the theoretical and empirical scenarios have recently seen new predictions 
based on Synthetic Horizontal Branch (SHB) models (Cassisi et al. 2004); Catelan, Pritzl, 
& Smith 2004) and new data (Butler 2003; Borissova et al. 2004; Storm 2004), we still 
lack a detailed NIR photometric analysis of a stellar system which includes a large sample 
of fundamental (F, RRab) and first-overtone (FO, RRc) RR Lyrae stars (only a dozen 
Galactic GCs host more than 50 RR Lyrae, sec, e.g., Clement ct al. 2001) and which 
spans a broad range in metal abundances. At present, the largest NIR data set includes 
74 RR Lyrae stars whose metal abundance, based on 23 objects, ranges from [Fe/H]pa —2 
to [Fe/H]f» —0.7, as observed by Borissova et al. (2004) in the inner regions of the Large 
Magellanic Cloud. In this context, the GC u; Gen is a perfect target, since it includes more 
than 186 RR Lyrae (86 RRab , 100 RRc ; Kaluzny et al. 2004) with a metal content ranging 

from [Fe/H] 2 to [Fe/H] 1 (131 objects, Rey et al 2000; 74 objects Sohima et al. 

2006). Moreover, the uj Cen distance has been estimated using different distance indicators 
such as the First Overtone Blue Edge (FOBE, Caputo ct al. 2002), the eclipsing binary 
OGLEGC-17 (Thompson et al. 2001; Kaluzny et al. 2002), the y-band metallicity relation 
of RR Lyrae stars (Rey ct al. 2000; Catelan 2005), and a dynamical analysis by van de Ven 
et al. (2006) based on proper motion and radial velocities of a large sample of individual 
stars. 

In this paper, we present new NIR photometry for a substantial number of the RR Lyrae 
stars in u Cen and we provide new distance estimates which we compare with determinations 
already present in the literature. 

2. Observations and data reduction 

The present NIR data set includes three different samples of J, Kg data either collected 
with the NIR camera SOFI available at the New Technology Telescope (NTT, ESQ, La Silla), 
or by Longmore et al. (1990), or available from the Two Micron All Sky Survey (2MASS) 
catalogue^. 

Near-Infrared J and Kg images of uj Cen were collected in two different runs, namely 
6,8 February 2001 (UT) and 3,4,25, and 26 February 2002 (UT). The pixel scale of SOFI is 
0.292" while the field of view is 4.94' x 4.94'. We observed three different fields: Field A is 



^See http://www.ipac.caltech.edu/2mass/releases/allsky/ 



-4- 



centered on the cluster center, field C at a distance ~ 10.5' NW away from the center, and 
field D at a distance ~ 10.7' SW away from the center (sec solid rectangles in Fig. 1). We 
collected 36 J- and 49 i^'g-band images of field A, 33 J- and 52 i^T^-band images of field C, 
and 12 J- and 18 Xg-band images of field D. For J images the exposure time was t = 3s, 
while for Kg images was t = 12 s. Owing to the crowded nature of the fields the observing 
time was equally divided between observations on uj Cen and sky observations. Together 
with these data we retrieved from the ESO archive a mosaic of nine J and nine Ks images 
of to Cen (SoUima et al. 2004) which cover an area of ~ 13' x 13' across the cluster center 
(see dashed rectangles in Fig. 1). Note that our Field A overlaps with the central pointing 
of this dataset. These images were collected on 13,14 January 2000 (UT), using the same 
experimental equipment. Seven out of the nine fields were observed once in J, Kg bands, 
while the other two fields were observed twice. The individual exposure times were 54 s and 
180 s. 

The pre-processing was performed using standard IRAF procedures. To improve the 
photometric accuracy we did not stack the individual images, but determined the (linearly 
variable) PSF for each frame based on more than 50 PSF stars uniformly distributed across 
each image. Photometry for the entire set of J and Kg images was then simultaneously 
performed with DAOPHOT/ALLFRAME (Stetson 1994) together with a few V- and I- 
band images (Dall'Ora et al. 2004; Del Principe et al. 2005). A more detailed analysis 
concerning the reduction strategy will be discussed in a forthcoming paper (Del Principe et 
al. 2006, in preparation). 

The absolute calibration in the 2MASS system (Cutri et al. 2003) was performed using 
a large set of 400 local stars selected from the 2MASS point-source catalogue. To constrain 
on a quantitative basis the intrinsic accuracy of both absolute and relative zero-points, we 
performed an independent calibration using a set of 4 standard stars (Persson et al. 1998) 
observed during the night of 24 February 2002 at air masses which bracket the observations of 
fields A, C, and D. The stars selected were S791-C, S867-V, S273-E, and S870-T. Fig. 2 shows 
the comparison between secondary standards in field A and local 2MASS stars. The observed 
distribution of magnitude differences has a standard deviation of K2M~Kp = — 0.012±0.076 
and J2M — Jp = 0.003 ± 0.080, which is a measure of the typical difference for an individual 
star. Note that the quoted differences agree well, within the errors, with the transformations 
provided by Carpenter (2001). As a further independent test, we also cahbrated the NIR 
data using standard stars collected during the other observing nights, and we found that the 
accuracy of the absolute zero-points is ~ 0.02 mag. 

We ended up with a photometric catalogue of ~ 115, 000 stars with hmiting magnitudes 
of J 20 and Kg ~ 19.0 mag. Fig. 3 shows the NIR CMD of the stars with a 'separation 
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index' , sep^ larger than 3.0 and an intrinsic photometric accuracy of uj^k < 0.08 mag. 

We identified 56 catalogued (OGLEII, Kaluzny et al. 2004) RR Lyrae stars (27 RRah , 
29 RRc ) in fields A, C, and D (see circles in Fig. 1) and 58 additional RR Lyrae (25 RRah , 33 
RRc ) variables (see diamonds in Fig. 1) in SoUima's fields. In order to increase the sample 
of RR Lyrae stars we also included data for the 29 RR Lyrae stars (19 RRab , 10 RRc ; see 
squares in Fig. 1) observed by Longmore et al. (1990). We have transformed their mean 
i^'-band data from the Anglo Australian Observatory NIR system to the 2MASS Kg system 
using the transformations provided by Carpenter (2001). Additionally, we cross-correlated 
the catalogues provided by OGLEII and by Clement et al. (2001) with the 2MASS catalogue. 
In this way we identified 37 more catalogued RR Lyrae stars (10 RRab , 27 RRc ; see triangles 
in Fig. 1), and therefore we ended up with a total sample of 180 RR Lyrae {SlRRab , 99RRc ) 
for which we have at least one J or X^-band measurement. 

To obtain the most accurate mean magnitudes for the RR Lyrae, we fit the individual 
Kg phase points measured by ALLFRAME with a template curve (Jones et al. 1996) and 
compute on the basis of this curve the intensity-averaged mean magnitude. The method 
requires an accurately known ephemeris of the epoch of maximum for each star as well as 
the amplitude in the optical V or B bands. These parameters are available from Kaluzny 
et al. (2004)^ for 167 of the variables in our sample. For each of the remaining 13 variables, 
all with single epoch 2MASS measurements, the 2MASS measurement was adopted as the 
best estimate of the mean magnitude. 

The typical accuracy of mean K — band magnitude of RR Lyrae stars either with a 
good coverage of the light curve or with single epoch measurements, epoch of maxima and 
luminosity amplitudes is better than 0.02 mag. In order to constrain the typical uncertainty 
for the remaining thirteen objects (2MASS sample) with only single epoch measurements, 
we estimated the difi^erence between the mean K — band magnitudes of the RR Lyrae stars 
for which we have a very good coverage of the light curves (22 RRab , 22 RRc ) and the 
mean K — band magnitudes only based on the 2MASS measurements. We found that the 
individual differences are smaller than 0.2 mag for RRab and 0.1 for RRc . 

The current RR Lyrae sample includes 16 Blazhko RR Lyrae stars (Kaluzny et al. 2004). 
Ten Blazhko RR Lyrae (eight in our sample and two in the Longmore sample) present a good 
coverage of the fight curves. For four objects (three in the SoUima sample and one in the 



^The 'separation index' estimates the photometric contamination due to crowding (Stetson et al. 2003). 
Current sep value corresponds to stars having less than 6% of their observed light contributed by known 
neighbours. 

^For more details see the home page of the Cluster AgeS Experiment, http://case.canik.edu.pl 
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2MASS sample) are available single epoch measurements and accurate epochs, and therefore 
the mean magnitude was determined using the template curve. For only two objects in 
the 2MASS sample we adopted the single epoch measurements, since recent epochs are 
not available. The Blazhko effect shows up as a modulation in the luminosity amplitude 
(Stothers 2006, and references therein) that might affect the accuracy of the mean K — hand 
magnitudes estimated with the template curve. In order to account for the uncertainty in the 
luminosity amplitude we estimated the mean K — hand magnitude by increasing/decreasing 
the amplitude by a factor of two. The difference in the mean K — hand magnitude is of 
the order of 0.02 mag. The error for the two objects with single epoch J and K — hand 
magnitudes is at most of the order of 0.15-0.20 mag. 

Unfortunately, no empirical template is available for the J-band light curves. Instead 
the mean J-band magnitude of RR Lyrae stars with well-sampled light curves (56 RR Lyrae 
stars located in fields A,C,D) was determined from a cubic spline fit to the hght curves. 
The fit was performed using spline under tension (Stellingwerf 1978, and references therein). 
This approach presents several advantages when compared with the classical Fourier fit, in 
particular for noisy light curves (Del Principe et al. 2006). For the other objects the J-band 
magnitude is based on a single epoch measurement. We found that the typical difference 
between the mean J — hand magnitudes of the RR Lyrae stars for which we have a very 
good coverage of the light curves (22 RRab , 22 RRc ) and the mean J — hand magnitudes 
only based on single epoch 2MASS measurements is systematically smaller than 0.3 mag for 
i?i?aband 0.15 mag for i?i?c stars. 

Open and filled circles in Fig. 3 show the location in the CMD of the entire sample of 
RR Lyrae stars. This is the largest sample of homogeneous NIR magnitudes for RR Lyrae 
stars in a stellar system ever collected. 

3. Distance determinations 

To determine the distance to u Cen on the basis of the observed mean dcrcddcncd Kq 
and Jo magnitudes we need to adopt a relation between observational quantities and the 
absolute magnitude. From the theoretical point of view, these relations have been recently 
investigated (e.g.. Bono et al. 2001, 2003; CPS04; C04) on the basis of Synthetic Horizontal- 
Branch (SHB) models with promising results in good agreement with empirical calibrations. 
These model grids can be parametrized in different ways, the simplest relation being the 
logP — Mk relation and the more complicated relations taking into account either the 
global metallicity Z (e.g. Bono et al. 2001, 2003) or the Lee (1990) Horizontal Branch type 
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thb^ (e.g., C04, and CPS04). 

Bono et al. (2003) showed that metaUicity has a significant impact on the i^-band 
magnitude of RR Lyrae stars. C04 showed that using the HB type as a parameter, at fixed 
metal abundance, also results in tight relations in both J and K. The latter approach is 
only valid for globular clusters where the HB type can be determined, and it is not directly 
applicable to field stars, uj Cen is a peculiar object as it is morphologically a globular cluster 
but the chemical composition and probably the ages of the stars differ significantly, making 
the population more similar to that of a dwarf galaxy. Consequently, the HB type is not 
well defined a priori, although the color-magnitude diagram presented in Fig. 3 shows that 
there are very few HB stars redder than the RR Lyrae stars. This feature suggests that all 
the sub-populations which contribute to the HB must have a collective HB type close to 1^. 

In order to determine accurate distance estimates for both F and FO pulsators, we 
computed a new set of SHB models with selected metal abundances (Z=0.0001, 0.0003, 
0.001, 0.004) and a helium content of y = 0.245. We did this by adopting the detailed set 
of HB models recently constructed by Pietrinferni et al. (2006) and the procedure already 
discussed by C04. The pulsation models discussed by Di Criscienzo et al. (2004) were used 
to fix the boundaries of the RR Lyrae instability strip and to evaluate the periods of the 
pulsators. On this basis, evolutionary and pulsation predictions were used to constrain the 
stellar distribution along the HB, and to compute for each individual simulation the HB- 
type Thb- The SHB models were constructed by accounting for the hysteresis effect (see, e.g.. 
Bono et al. 2005) and by neglecting the pre-ZAHB evolutionary phases (see, e.g., Piersanti, 
Tornambe', Castellani 2004). 

Theoretical predictions were transformed into the Bessell & Brett (1988) NIR photo- 
metric system using the bolometric corrections and color-temperature relations provided by 
Pietrinferni et al. (2004). They were then transformed into the 2MASS system using the 
transformations provided by Carpenter (2001) and eventually, for each given SHB, the pre- 
dicted Period-Luminosity relations, as given in the form Mj^k — o,j^K + bj,K ^og P + cj^kTrb^, 
were derived. The coefficients and the relative errors of the least-squares solutions are fisted 
in Table 1. The zero-points and the coefficients of the above PLZ relations further support 
the use of independent relations for F and FO pulsators. As a matter of fact, the slopes and 



''thb is defined as the ratio t^^b = (B — R) / {B + V + R), where B, V, R refers to the number of blue, 
variable, and red HB stars, respectively, thus resulting in the = — 1 for HBs with only red stars and thb = 1 
for HBs with only blue stars. 

^Castellani ct al. (2006, in preparation) find thb = 0.94 on average for lu Cen. 

^Note that in the following we will use -fsT-band notation for (2MASS)-band. 
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the zero-points between the two different groups differ at least at ten a level for both J and 

K-band. The referee suggested to check the dependence of current F and FO PL relations 
on the predicted effective temperature of F blue edge and FO red edge. We performed a few 
numerical experiments artificially increasing/decreasing by 150 K the effective temperature 
of FO red edge and of F blue edge. We found that the change in the effective temperature of 
the F blue edge affects the coefficients of J and /T— band FO PL relations by ±0.02 — 0.03 
(a) and by ±0.02 — 0.05 (6), respectively. The coefficients of J and band F PL relations 
change by ±0.02 (a) and by ±0.02 — 0.04 (6), respectively. The change in the effective tem- 
perature of the FO red edge affects the a, h coefficients of F and FO relations on average by 
±0.01. The impact of these changes on the coefficient of the HB type (c) is marginal. Note 
that the adopted uncertainty on the location of FO red edge and F blue edge is supported 
by the comparison with RR Lyrae stars in GCs (Di Criscienzo et al. 2004). These findings 
further support the approach to derive different PL relations for F and FO pulsators. It is 
also noteworthy, that the use of fundamentalized FO periods causes an increase of a factor 
of two in the uncertainty of their zero points and slopes. 

As a whole, these relations agree quite well with those already presented by C04 (see 
their Tables 6 and 8). In particular, i) the slopes of the Mj/^ — logP relations are, at 
fixed metal content, independent of the HB-type, while they become slightly steeper with 
increasing Z. ii) The zero-points of the Mj^k — log-P relations become, at fixed metal con- 
tent, brighter when moving from blue to red HB populations. This dependence significantly 
decreases with increasing Z. in) The zero-points of the Mj^k — log-P relations, indepen- 
dently of the HB-type, become fainter with increasing Z. An even more striking agreement 
is found with the relations computed by CPS04. For the relevant mean metallicity model, 
Z = 0.001, and HB type, thb = 0.934, they find for their combined sample of fundamen- 
tal mode and first-overtone pulsators = —2.388 log — 1.133 whereas we find for the 
fundamental mode pulsators = —2.37 log Pp — 1.13 and for the first-ovcrtonc pulsators 
Mk = —2.41 log Ppo ~ 1.20. The referee noted that, according to CPS04, the coefficients a 
and b of the PL relations present some dependence on both the HB type and the metal abun- 
dance. Unfortunately, a global comparison is hampered by the different approach adopted to 
derive the PL relations, i.e. multilinear least-squares vs. third order polynomials (see their 
equations. 1 and 2 and their Table 9). Moreover, the metallicity of current SHB models 
range from Z=0.0001 to Z=0.004, while those from CPS04 range from Z=0.0005 to Z=0.006. 

The individual reddening values of RRab stars have been empirically estimated using a 
Period-Color(y — ir)-Amplitude(Ay) relation (Piersimoni et al. 2002, 2006, in preparation), 
while for RR^ stars and RR Lyrae affected by amplitude modulation we used a reddening 
map based on u,v,b,y Stromgren photometry (Calamida et al. 2005). Accurate mean V- 
band magnitudes for RR Lyrae in our sample have been measured either by Kaluzny et al. 
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(2004)''' or by Butler et al. (1978) for the outermost ones. For the remaining objects, we 
adopted a mean reddening of E{B — V) = 0.11 ± 0.02. Selective absorptions have been 
estimated using the reddening law from CardcIIi ct al. (1989), that is Av/E[B — V) — 3.1, 
Aj/E{B -V) = 0.868, and Ak/E{B - V) = 0.341. 

Data plotted in Fig. 4 show the fit between the dereddened K magnitudes and predicted 
PL relations at fixed metal abundance and for thb = 0.94. In this figure, the asterisks mark 
the ovcrluminous RRah stars which are either blended with a fainter close companion (V32, 
VI 18, V139, and V150) or have a controversial mode identification (V84), while the open 
triangle refers to V182 which, together with three other (V168, V181, V183) undcrluminous 
variables with Kq >14 mag, presents peculiar J ~ K colors. By excluding these variables, we 
estimate a true distance of (/Xq) = 13.80 ±0.10 mag and 13.70 ±0.10 mag with Z=0.0001 and 
Z=0.001, respectively, clearly indicating that individual metallicity estimates are required 
to improve the distance determination. 

Conveniently enough, individual metal abundances on the scale of Zinn & West (1984), 
as based on the hk index, for a substantial fraction (130) of the RR Lyrae plotted in Fig. 
4 have been measured by Rey et al. (2000). In a recent investigation based on medium 
resolution spectroscopy, SoUima et al. (2006) found, for 74 RR Lyrae in u Cen , no sys- 
tematic difference with the metal abundances based on the hk index. To properly use the 
theoretical predictions, the observed [Fe/H] values have been transformed into "global" Z 
values by taking into account possible a-element enhancements'^. Following to Gratton, Sne- 
dcn, & Carctta (2004), who found [Q;/Fe]~ 0.45 as a plausible upper limit to the a-element 
enhancement, we adopt / = lO^"/-^^' = 3. 

We are now in a position to derive distance moduli of individual RR Lyrae stars with 
metallicity measurements from Rey et al. (2000). In panel (d) in Fig. 5 we show the 
resulting values as a function of metallicity, leading to {^q) = 13.70 ±0.06 mag. To illustrate 
the sensitivity of the method to possible errors in the assumed HB type or a-enhancement, 
we have also computed the distance moduli for trb = and/or f — 1 and plotted the results 
in Fig. 5. We note that the effect on the distance moduh is rather small considering that the 
estimated error on the HB type, cr^ = 0.2 (see star counts of HB stars in NGC 2808 listed in 



^For more details see the home page of the Cluster AgeS Experiment, http://case.camk.edu.pl 

^Following Salaris, ChiefR & Straniero (1993) the global metallicity was estimated as 
log Z= [Fe/H] +logZ0 + log(0.362 + 0.638/), where Zo=0.012 (Asplund et al. 2004) and / is the en- 
hancement factor of a-elements with respect to iron. Note that the change in the distance modulus 
when moving from scaled solar to a enhanced predictions also includes current uncertainties in the solar 
metallicity. 
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Table 2 of Castellani et al. 2006), and on the ct-enhancement factor, cr/ = 0.5, are typically 
smaller than the values adopted here. Of particular interest is the fact that the effect of an 
error in the assumed HB-typc is to introduce a slope in these diagrams (see panels a and 
b). According to the relations given in Table 1, this trend can be explained as an empirical 
evidence that the more metal-poor RR Lyrae are connected with bluer HB populations when 
compared with the more metal-rich ones. 

The i^'-based distance determinations are fully supported by the J-band measurements. 
The estimated uncertainty on the mean J-magnitudes depends on the luminosity amplitude, 
{aj{RRc ) ~ 0.1; (jj{RRab ) ~ 0.2), and it is less accurate than for the K-hand measurements 
{ax ~ 0.015). It is worth mentioning that we still lack empirical template curves for the 
J — band, and therefore accurate mean J — band magnitudes are only available for 56 out of 
the 180 RR Lyrae in our sample. As shown in Fig. 6, the true distance moduli based on the 
relations given in Table 1 agree with the values based on X-magnitudes, but with a larger 
error. Although, current mean J — band magnitudes present this drawback, this NIR band 
appears very promising to estimate distances, since it presents several advantages over the 
K — band and requires much less telescope time to reach similar photometric accuracy. 

We can also determine the distance to the individual RR Lyrae stars using the PLZ/^ relations 

of Bono ct al. (2003). By using both field and cluster RR Lyrae stars they found Mk — 
—0.77 — 2.101 logPp -|- 0.231 [Fe/H] where the fundamcntalized logarithmic period was com- 
puted as log Pp = log Ppo + 0.127. Combining this relation with the sub-sample of RR Lyrae 
for which are available metallicity estimates from Rey et al. (2000), we find /xq = 13.77±0.07 
mag based on a total of 113 RRabSmd i?i?c stars. 

In addition to the quoted intrinsic errors we estimate the external, systematic error 
to be (Jsys = 0.06 mag. This is the result of summing in quadrature the contributions 
from the photometric zero point {a = 0.02 mag), the uncertainty on the absorption 
{a = 0.01 mag, based on uncertainties of 0.1 on Rk and 0.02 mag on E{B — V)). The 
uncertainty of the reddening is based on the recent Stromgren photometric investigation by 
Calamida et al. (2006), while the uncertainty on the selective absorption coefficient is due to 
current uncertainties in the adopted extinction law (Rieke Sz Lebofsky 1985; CardeUi et al. 
1989, and references therein). The absolute zero point of the metallicity scale is a significant 
source of uncertainty as discussed by Dall'Ora et al. (2004) due to systematic differences 
between metallicity scales. They estimate a = 0.25 dex on the zero point which results in 
the dominating contribution to the systematic error of cr = 0.25 x 0.231 = 0.06 mag. 

As the LV Cen RR Lyrae stars exhibits a broad range of metallicities and all belong to 
the same cluster, we can empirically constrain the metallicity effect without having to rely 
on uncertain distance estimates. We have made linear three-parameter fits of the relation 



in Eq.l for RRabSind RRc stars separately, using only dereddened apparent magnitudes, and 
using only the variables with individual hk abundance evaluations. We find 



Ko{F) 


= 12.70(±0.06) 


- 2.71(±0.12) logPp + 0.12(±0.04)[Fe/H] 


{N = 


54) 


(1) 


Ko{FO) 


= 12.26(±0.05) 


- 2.49(±0.15) log Pfo + 0.04(±0.04)[Fe/H] 


{N = 


48) 


(2) 


MF) 


= 13.09(±0.12) 


- 2.37(±0.25) logPi. + 0.15(±0.08)[Fe/H] 


(N^ 


54) 


(3) 


UFO) 


= 12.58(±0.08) 


- 2.09(±0.21) log Pfo - 0.01 (±0.05) [Fe/H] 


{N^ 


48) 


(4) 



where the symbols have their usual meaning. In Fig. 7 we have plotted the observed 
magnitudes against metallicity after having removed the period dependency. The metallicity 
effect is weak for the fundamental pulsators and for the first-overtone pulsators the metallicity 
effect is consistent with only a marginal dependence. In both cases the relations are in 
reasonable agreement with the predictions from the SHB models. Current K — band PL 
relations present two substantial differences when compared with the PLZ relation derived 
by Bono et al. (2003): a) we derived independent PL relations for F and FO pulsators, while 
Bono et al. only derived semi-empirical "average" PL relations for F and fundamentafized FO 
pulsators; b) current PL relations only apply to cluster variables for which we can estimate 
the HB type, while the PLZ relation derived by Bono et al. (2003) appfies to individual RR 
Lyrae stars. 



4. Conclusions and final remarks 



In Table 2 we have listed recent empirical distance estimates to u Cen where we have 
corrected all the moduli for the same reddening and employed the same reddening law as 
we have used for the RR Lyrae stars. We note that our new near-IR results are in excellent 
agreement with the largely geometric distance to the eclipsing binary OGLEGC-17 from 
Thompson et al. (2001) and Kaluzny et al. (2002), and from the First Overtone Blue 
Edge method based on pulsation predictions from Caputo et al. (2002). Current distance 
determinations based on J and K — band measurements and on NIR PL(rHB) relations 
are also in very good agreement with similar estimates based on NIR relations derived by 
CPS04. In particular, for thb = 0.94 and / = 3 we found 13.70 ± 0.06 vs 13.72 ± 0.06 
{K — band) and 13.71 ± 0.10 vs 13.76 ± 0.10 (J — band). Note that distance estimates based 
on CPS04 relations were determined using the same solar metal abundance and we did not 
include a few very metal-poor RR Lyrae, since CPS04 prediction range from Z=0.0005 to 
Z=0.006. This finding is also supported by most recent My— [Fe/H] relations. By adopting 
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the mean V-magnitudes for RR Lyrae stars provided by Kaluzny et al. (2004) and by Butler 
et al. (1978) together with the calibration given by Bono et al. (2003) wc obtain for RR 
Lyrae with individual metal abundances a distance value of (m — M)o = 13.72 ih 0.11. We 
obtain a similar distance modulus (13.62 ± 0.11) using the same magnitudes and the recent 
calibration of the My— [Fc/H] relation provided by Catelan (2005). The slope of such a 
relation was estimated by Cacciari & Clementini (2003), while the zero-point is based on the 
trigonometric parallax of RR Lyr itself (Benedict et al. 2002; Bono et al. 2002). Note that 
the systematic error for distance estimates based on visual magnitudes is cr = 0.08, due to 
the significant contribution from the absorption correction. 

On the other hand, we deviate by more than three a from the distance evaluation based 
on High Amplitude S Scuti (HADS) stars (/xq = 14.05 ±0.02, McNamara 2000) and from the 
modulus of Ho — 13.41 ± 0.13 from van de Ven et al. (2006) based on cluster star dynamics. 
It is worth noting that current uncertainties on metal abundance and on NIR photometry 
can not account for such a discrepancy. 

The NIR PL relations for RR Lyrae stars provide an accurate approach for determining 
absolute distances to globular clusters and other stellar systems. We find for u; Cen a dis- 
tance modulus of 13.70 ± 0.06 ± 0.06 where the error estimates are random and systematic 
errors, respectively. We note that for populous systems like globulars and dwarf galaxies the 
dominating source of error are the systematic sources. The accuracy of NIR photometry in 
crowded regions significantly improved during the last few years and the sample sizes will 
certainly benefit from the new large field of view NIR cameras available at 4m class tele- 
scopes (van der Bliek et al. 2004). New J— band light curve templates for both fundamental 
and first overtone RR Lyrae are also urgently required to improve the accuracy of their mean 
magnitudes in this band. 
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Table 1. Predicted near-infrared Period-Luminosity relations. The pulsation modes are 
indicated as F for fundamental and FO for first-overtone pulsators. 
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Table 2. Distance determinations to u Cen. 



Method 


(m - M)o'^ 


colheadNotes 


PL-THB {K) 


13.70 ±0.06 


b 


PL-THB {J) 


13.71 ±0.10 


b 


PLZx 


13.77 ±0.07 


c 


PL-THB (^) 


13.72 ± 0.06 


d 


PL-THB iJ) 


13.76 ±0.10 


d 


Mv - [Fe/H] 


13.72 ±0.11 


e 


My - [Fe/H] 


13.62 ±0.11 


Catelan et al. (2005)^ 


Eclipsing binary 


13.71 ±0.11 


Thompson et al. (2001) 


Eclipsing binary 


13.75 ± 0.04 


Kaluzny et al. (2002) 


FOBE 


13.74 ±0.11 


Caputo et al. (2002) 


High ampl. 5 Set 


14.05 ±0.02 


MeNamara (2000) 


Cluster dynamics 


13.41 ±0.13 


van de Ven et al. (2006) 



^Using the same mean reddening (E(B — V) = 0.11 ± 0.02) 
and reddening law (Cardelli et al. 1989) adopted in this paper. 
The error estimates do not include the contribution of systematic 
errors (crj^g = 0.06 mag, crYys = 0.08 mag). 

''Based on the PL-thb relations presented here for thb = 
0.94, an a-elemeut enhancement factor / = 3, and the new solaj: 

metal abundance. 

"^Based on the semi-empirical K—hwd PLZ relation from 
Bono et al. (2003). 

''Based on the PL-thb relations from Catelan et al. (2004) 
and Thb — 0-934, an a-element enhancement factor / = 3, and 
the new solar metal abundance. 

''Based on the Mv-[Fe/H] relation from Bono et al. (2003) 
and mean V-magnitudes from Kaluzny et al. (2004) and Butler 
et al. (1978). 

'Based on the M^ — [Fe/H] relation from Catelan (2005) and 
mean V-magnitudes from Kaluzny et al. (2004) and Butler et 
al. (1978). 
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Fig. 1. — Distribution of RR Lyrae stars in the field of oj Cen. Open and filled symbols show 
first-overtone and fundamental mode RR Lyrae. Solid and dashed rectangles show the fields 
observed with SOFI@NTT. See text for more details. 
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Fig. 2. — Difference in K (top) and J (bottom) magnitudes of Field A based on two different 
absolute zero-point calibrations. Jpersson, Kpf^rsson secondary standards have been calibrated 
using four standard stars measured by Persson et al. (1998), while J2MASS, K2MASS sec- 
ondary standards have been calibrated using pa 400 local 2MASS standard stars. 
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Fig. 3. — NIR Color-Magnitude Diagram of u Cen. The stars plotted in this figure are 
located within ^ 13' x 13' across the cluster center. Open and filled circles mark the sample 
of first-overtone and fundamental RR Lyrae stars for which we estimated J and K-hand 
magnitudes. Error bars display intrinsic photometric errors. 
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Fig. 4. — Dereddened K magnitudes of RR Lyrae stars in u Cen versus period. Dashed 
and solid lines display the predicted Period-Luminosity relations for F (filled circles) and FO 
(open circles) pulsators, respectively. The adopted metal content and HB-type are labeled. 
The error bars show intrinsic errors in the measured magnitudes. 
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Fig. 5. — Individual true distance moduli versus [Fe/H] for RR Lyrae in cu Cen. The panels 
show distance evaluations for scaled solar (/=!) and a-enhanced (/=3) chemical composi- 
tions, as well as with HB-type thb = and 0.94. The error bars in the bottom right panel 
display individual errors affecting K magnitudes and [Fe/H] abundances. The error bars 
in the other panels have been omitted to make more clear the distribution of F and FO 
variables. 



-23- 



14.4 
14.2 



a, 



14 - 



:l 13.8 

13.6 

13.4 
14.4 

14.2 
1 14 

N 

13.8 
13.6 
13.4 



I I I I I I I I I I I I I I I I I I I I I I I I 

a) f=l, T„B=0 
_ </io2> = 13.81±0.12 mag J 


























-++- 


-++++- 


-++++- 


-++++- 


-++++- 


+ 




1 1 1 1 


-++++- 


-++++- 


III 



b) f=3, T„3=0 

<At(,2> = 13.74±0.11 mag 




I I I I I I I I I I I I I I I I I I I I I 



c) f=l, T„B=0.94 
_ </i„2>=13.'''6±0.11 mag 




d) f=3, T„B=0.94 
<At(i2>=13.71±0.10 mag 




-2.5 -2 -1.5 -1 -0.5 
[Fe/H] 



-2.5 -2 -1.5 -1 -0.5 
[Fe/H] 



Fig. 6. — Same as in Fig. 5, but for the mean J-magnitudes. 



-24- 




_l I I I I I I I I I L_ 

-2 -1.5 -1 

[Fe/H] 



Fig. 7. — The mean Kq magnitude corrected for the empirically determined period effect as 
a function of metallicity for fundamental {RRab ) and first-overtone {RRc ) pulsators respec- 
tively. The best fit hues have been overplotted. 



